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We generalize the Landau-Zener theory of coherent tunneling transitions by taking thermal relax- 
ation into account. The evaluation of a generalized master equation containing a dynamic tunneling 
rate that includes the interaction between the relevant system and its environment leads to an inco- 
herent Zener transition probability with an exponent that is twice as large as the one of the coherent 
Zener probability in the limit T ^ 0. We apply our results to molecular clusters, in particular to 
recent measurements of the tunneling transition of spins in Peg crystals performed by Wernsdorfer 
and Sessoh [Science 284, 133 (1999)]. 
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The adiabatic transition in a two-level system 
{|m) , \m')} with energy level crossing is described by the 
Landau-Zener transition probabilityfil 
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where (e™' ~ ^m) is the energy of the state |m) (|m')), 
and E„im' is given by the coupling between these states. 
Equation (|l|) is obtained direcily from the Schrodinger 
equation with the Hamiltoniann 
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If the total Hamiltonian Titot comprising Ti. forms a po- 
tential barrier — as we shall assume from now on — 
Emm' corresponds to the tunnel splitting energy, Pcoh to 
the coherent Zener tunneling probability, and the eigen- 
states of Ti. are delocalized as long as em — £m' ^ Emm' ■ 
As there is a large amount of potential barriers in physical 
systems, Eq. (|l]) has becomejaii important tool for study- 
ing tunneling transitions .Ba'D'El It must be nated that all 
quantum systems to which the Zener modeld is applica- 
ble can be described by 'pure states and their coherent 
time evolution. It is the aim of the present work to gen- 
eralize the Zener theory in the sense that we take also 
the incoherent evolution of mixed states into account (see 
also Refs. ^ and ^, 0, ||, ^, |l^ for a comparison). In or- 
der to provide a clear description of our generally valid 
theory, we give the derivation of the incoherent Zener 
tunneling probability Pine (see Fig. ^) in the framework 
of spin tunneling in molecular magnets, which has be- 
come a highly attractive research field in the past few 
years since several experiments revealed mesoscopically 
observable quantum phenomena in molecular clusters, 
such as Mni2-acetate (Mni2) (Refs. 
[l6| ) and Fes-triazacyclononane (Fes)) 
ticular, we will show that our theory presented in this 
work is in good agreement with recent measurements of 
Pine as a function of the extqs^ri transversal field for 
various temperatures in Fcg 




FIG. 1. Energy level crossing diagram for incoherent Zener 
transitions. Dotted lines: transitions due to interaction with 
environment, leading to a linewidth '^rnm' ■ Variables are ex- 
plained in the text. 

We proceed now from the assumption that the range 
over which £mm'{t) = £m ~ £m' is swept, defined by 
the boundaries e<„j, := miiitiemm'} and e>,„, := 
maxtjemm'}, is much larger than Emm' and the deco- 
herence rate h'jmm' (see below and Fig. ^). In addition, 
we restrict the evolution of our system to times t that are 
much longer than the decoherence time Td = 1 /^mm' ■ In 
this case, tunneling transitions between pairs of degener- 
ate excited states are incoherent. This tunneling is only 
observable if the temperature T is kept well below the 
activation energy of the potential barrier. Accordingly, 
we are interested only in times t that are larger than 
the relaxation times of the excited states. Thus, we can 
apply our formalism presented in Ref. ^ which treats 
incoherent tunneling between pairwise degenerate states 
within a single spin system, the decoherence of which is 
due to the interaction with its environment. We showed 
in Ref. ^ that one can reduce the generalized master 
equation comprising off-diagonal elements of the density 
matrix p to a complete master equation that consists only 
of the diagonal elements, i.e., 

Pm = -WmPm + ^ Wm7iPn + F" {pm' - Pm) , (3) 
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is the incoherent tunnehng rate from \m) to |rn'), which, 
in contrast to Ref. is assumed now to be time de- 
pendent (see below for range of validity). In Eqs. (H) 
and we have made use of the abbreviations 7m,m' = 

{Wm + Wra')/2 and Wm = Y.n ^nm, whcrC Wnm dcuOteS 

the approximately time- independent transition rate from 
|m) ±D \n), which can be obtained via Fermi's golden 
ruleBj For straightforward calculation of Emm' it is u^ 
ful to note that our generalized tunnel splitting formulate 
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can be represented by the graphs shown in Fig. ||. Vmi,mj 
denote off-diagonal matrix elements of the total Hamil- 
tonian Titot ■ 
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FIG. 2. Tunnel splitting energy Emm'- The open circles o 
correspond to the states \mi), i = 1, . . . , N , the solid circles • 
to |m), I'm'), and the lines to the matrix elements Vm.i,mj ■ 

First we solve Eq. in the unbiased case — cor- 
responding to n = (see below) — where the ground 
states |s), |— s) and the excited states \m), |— m), m G 
[[s] — s -I- 1, s — 1] of our spin system with spin s are pair- 
wise degenerate. In addition, we assume that the excited 
states are already in their stationary state, i.e., pm = 
Vto 7^ s, — s. This implies that Smm'it) in Eq. (^) must 
be changed within a time that is much smaller (adiabatic 
approximation) or much greater (sudden approximation) 
than the relaxation times of the excited states, which are 
of the order oil /Wm- Proceeding as in Ref. (Sec. V A) 
we obtain from Eq. (H) 
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where we have defined the quantity Ap(t) ^ pa — P-s, 
which satisfies the initial condition Ap(i = to) — li 
and thus Pinc(* = ^o) = 0. In distinction to Eq. (|^), 
we call Pine the incoherent Zener transition probabil- 
ity. The total time-dependent relaxation rate is given by 
Ftot = 2 [rj** -I- Fth], where the thermal rate Fth, which 
determines the incoherent relaxation via the excited, 
states, is evaluated by means of relaxation diagrams,E3 
such as shown in Fig. |^. For example, if we allow 
only for thermal transitions with Am = 1, we obtain 
Tth = /s-i ° /s-2 o • ■ • o /[s]/s(0). This continued fraction 
is recursively defined by 
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FIG. 3. Unbiased (n = 0) relaxation diagram of a 
spin system with spin S and symmetric anisotropy barrier 
7ia ~ —ASl- The solid lines correspond to thermal transi- 
tions with Am G [1,2s], double arrows indicating that there 
is more than one incoming rate, and the dashed lines represent 
tunneling transitions. 

Assuming linear time deuendence, i.e., Smm'{t) = 
a™ in the transition region jd and with je^j',^/ 1 ^ fi-'lmm' 
we obtain from Eq. m) 
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where we have set to — —t. In the low-temperature limit 
T — > the excited states are not populated anymore and 
thus Fth, which consists of intermediate rates that are 
weighted by Boltzmann factors 6m,E3 vanishes. Conse- 
quently, Eq. (^ simplifies to 



Ap = exp ■ 
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where the second expression is more general and can also 
be obtained directly from Eq. (0) by reducing 
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We note that the exponent in Eq. (g) differs by a factor 
of 2 from the Zener exponent in Eq. (|l|). This is not 
surprising since Ftot is the relaxation rate of Ap, where 
both ps and p-s are changed in time by the same amount, 
and not an escape rate like in the cases of coherent Zener 
transition and a decay, where only the population of the 
inital state is changed in time. Note that Eq. (||) implies 
^inc = 1 for |es._s(0)| (adiabatic limit) and Pine — 
for |es,-s(0)| — > oo (sudden limit). 

Instead of the linear time dependence es,-s — ctj'^t, 
one can consider oscillations of the form es,-s = asinujt. 
Neglecting Fth in Eq. (0) in the limit T ^ we get 
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where we have set rj = ^ a? -\- tfl^fl _g. Integrating 

Ftot from to = — 7r/2w to t = 7r/2w we obtain Ap — 
exp |— 7ri?^ _g/?i(jj?7}. In comparison to Eq. (||) we get 
here the extra factor 77, which provides an experimen- 
tally exploitable dependence on 7s, _s, provided that 
a < h-js-s- 

If we apply a bias to our system in such a way that our 
states are tuned to other resonances (see Fig. e.g., n = 
1,2,... (see below), the total relaxation rate is changed 



to rf-« « 2/ [l/(r7^+i +rj^-) + l/W^_,,_,+iJ, where 



the approximation comes from the fact that there can be 
other nonvanishing thermal rates W-s,n with initial state 
|n), n 7^ — s + 1, and final state |— s). determines 
the relaxation through the states that have higher energy 
than Is). 
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FIG. 4. Example of a biased relaxation diagram of a spin 
system with spin Sz and asymmetric barrier 7ia + Ti-z for 
the resonance condition n — 1. The diagram is explained in 
Fig. |. 



his paper, we apply our theory to 
'^^ which measured quantum oscil- 



In the second par 
recent experiments, 
lations of the tunnel splitting Emm'{Hx) in Fes as a func- 
tion of an externally applied transversal magnetic field 
Hx ■ In a coherent spin-state path integral approach these 
quantum oscillations and their associated spin parity ef- 
fects can be viewed as a result of interfering Berry phas 



carried by spin tunneling paths of opposite wi 
which are modified in the presence of a field 
ever, for a quantitative analysis of the Feg dat 
operator formalism presented here proves to be more use- 
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ful than the path integral approach.— ___ _ 
In accordance with earlier worHaliaHBH El 
single-spin Hamiltonian H = Ha + Ht + Hz 
describes sufficiently well the behavior of the giant spin 
S with s = 10 of a Fes cluster. It turns out ttwir-aur 
theory is in optimal agreement with experimentsE2IEil if 
we choose the easy-axis Hamiltonian to be Tirs-,= —AS'^, 
with anisotropy constant A/ks = 0.275 kJ13 and the 
in-plane Hamiltonian to be 
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with the anisotropy constants -62/^5 = 0.046 Kji3'Ea 
B^/kB = -6.0 X 10"^ K, Be/ks = 2.0 x 10"^ K, and 



Bg/ks = 2.0 X 10~^^ K. Besides the B4 term introduced 
in Refs. |2^ and |2| to obtain the desired period, we find 
that the Bq term is necessary to achieve the desired tun- 
neling amplitude in Fig. ^, while the Bs term is respon- 
sible for the minimum at ~ 1.4 T. The Zeeman cou- 
pling g/isS • H has been divided into a longitudinal part 
Tiz = gfJ-sHzSz and a transversal part, being the second 
term in Eq. (11), where H is the magnitude of the exter- 
nal magnetic field H, and z9 and ip define the spherical an- 
gles. According to Eq. (||), Ht induces tunneling between 
pairwise degenerate Sz eigenstates |m), —s < m < s, of 
Hn + Hz, with eigenvalues e„i. The resonance condition 
for such degeneracies, i.e., Em — £m', leads to the res- 
onance field JJ™"* = nA/g^B, n = m + m' . As can 
be seen in Fig. ^ our analytic formula (|^) for the result- 
ing tunnel splitting Emm' {Hx) is in reasonable agreement 
with the exact diagonalization of Ha, + Ht + Hz , which 
provides a good fit of the data in Ref. EOl 
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FIG. 5. The tunnel splitting energy E\o^-io{Hx) in Fes, 
exhibiting Berry phase oscillations, was calculated by exact 
diagonalization (solid line) and by using the approximate an- 
alytic formula (Bl) (dashed line). The angles are i9 — 90°, 
Lp = 4° . The period is in excellent agreement with data 
(Refs. 1^ and ^ if we set = 1.9 (Ref. The tun- 

nel splittings for — read: _Eio,-io ~ 9.0 x 10^* K, 
£9,-9 = 6.5 X 10"® K, and Es.-s. = 2.1 x 10"'' K. 

In order to account for thermal transitions between 
the \m) states, we include the most general spin-phonon 
couplingc3 which is allowed ia leading order by the D2 
symmetry of the Fes crystal,cZl i.e., 
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Hsp — giexxSx + g2>^vy-y 

+ gA^xz {Sx,Sz} + gh^yz {Sy, Sz} 
+ geojxy {Sx, Sy} + gjLUxz {Sx, Sz} 

+ gs(^yz {Sy,Sz}) , 
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where gi, i = 1,...,8, are the spin-phonon coupling 
constantSj— which we assume to be approximately equal, 
gi « ffo-l23 We knowl23 that within the spin system 
the first- and second-order thermal transition rates 
Wm±i,m, Wm±2,m are the strongest ones, which are eval- 
uated by Fermi's golden rule to be (r = 1, 2) 
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where s±i = (s =F to)(s ± m + l)(2m ± 1)^, s±2 = (s =F 
m){s ±m+ l)(s ^m — l)(s ± m + 2), and gi = 48, (72 ~ 
32. The^mass density p for Fcs is given by 1.92 x lO'^ 
kg/m'^,Lj the sound velocity c by 1400, m/s (yielding a 
Debye temperature of 0d — 33 K).Ea The-, incoherent 
Zener probability Pine in Fig. |^ fits the datacil well if one 
adjusts the couphng constant to go = 2.3 K. 
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FIG. 6. Unbiased (n — 0) relaxation diagram for Feg. Solid 
(dashed) lines: thermal (tunneling) transitions. 

For the temperature range 0.05 K< T < 0.7 K we 
achipse good agreement between our theory and the 
dataEJ if we take the states |±10), |±9}, and |±8) into 
account. In particular, the path leading through |±8) 
gives a non-negligible contribution for T > 0.6 K. Solv- 
ing the relaxation diagram shown in Fig. ^ we obtain from 
Eq. (^ the following results for Feg in the case n = 0: 
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FIG. 7. Zener transition probability Pinc{Hx) for temper- 
atures T =0.7 K, 0.65 K, 0.6 K, 0.55 K, 0.5 K, 0.45 K, and 
0.05 K. By choosing B4, = -6.9 x 10"'^ K for this plot our 
fit agrees well with data (Ref. except for the minima at 
Hx ~ ±0.2 T, which are too narrow (see text). Note that 
Pine is equal to 2P in Ref. ^ The tunnel splittings for this 
figure read: Pio,-io = 1-3 x lO"'^ K, Pg-g = 8.8 x 10"® K, 



and Pg 



= 2.7 X 10"* K for = 0. 
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where we have used the approximation 75 
and I £<'>,! > 



n ~ Wio.n 

P«,-n,7n.-ri- Pine = 1 " Ap, which 

is plotted in Eig. 0, is in good agreement with the 
measurements,Ej except for the most narrow minima at 
Hx w ±0.2 T. However, the experimental uncertainty of 
these minima is very largp-siace they depend strongly on 
the initial magnetizationEd'Ej In order to account for the 
increase of Pine at T — 0.6, 0.65, 0.7 K for higher fields 
Hx, we had to correct the energy levels Em occurring in 
Eq. (|l|) by first -order perturbations in H^;. Below 0.4 
K, Pine is T independent, which agrees well with Ref. |2^. 

In conclusion, our theory, which is based only on 
thermal-assisted tunneling and neglects dipolar and 
hyperfincEj fields, agrees well with the recent measure- 
ments in Refs. |o|and 21 and also leads to the prediction 
of the anisotropies P4, Bq, Pg, and the spin-phonon cou- 
pling constant go in Feg. 

Detailed calculations of the biased cases (n = 1,2,...) 
will be published elsewhere. 
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' Hyperfine effects in Fes become noticeable for T < 1.5 
K, and they lead to measurable corrections to Pine, which 
are small in the sense that the overall shape of the curve 
finc(T) is well described by spin-phonon effects (Ref. Uw- 
Also we note that the value for i5io,-io reported for ^ Fes 
and for ^^Feg differ by less than 20% (Ref. ^ , whereas the 
value reported initially for '''Feg (Ref. and the one re- 
ported subsequently for "^Fes (Ref. ^l]) differ by more than 
25%. Given this experimental uncertainty further measure- 
ments would be desirable. 
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